Biphoton ququarts as either pure or mixed states, features and reconstruction from 

coincidence measurements 
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Features of biphoton polarization-frequency ququarts are considered. Their wave functions are 
defined as functions of both polarization and frequency variables of photons with the symmetry 
obligatory for two-boson states taken into account. In experiments, biphoton ququarts can display 
different features in dependence on whether experiments involve purely polarization or (alterna- 
tively) polarization-frequency measurements. If in experiments one uses only polarization measure- 
ments, the originally pure states of ququarts can be seen as mixed biphoton polarization states. 
Features of such states are described and discussed in details. Schemes of coincidence measurements 
for reconstruction of the ququart's parameters are suggested and described. 

PACS numbers: 03.67.Bg, 03.67.Mn, 42.65.Lm 



I. INTRODUCTION 



Biphoton polarization-frequency ququarts can be pro- 
duced in processes of collinear Spontaneous Parametric 
Down- Conversion (SPDC) non-degenerate with respect 
to frequencies of photons [1]. In such states photons have 
two degrees of freedom: polarization and frequency. In 
terms of photon polarization and frequency variables, a 
and w, each of them can take independently one of two 
values: a = HoiV (horizontal or vertical polarization) 
and u) — u>h or (high or low frequencies) Q , @ • For 
experimental investigation of such states one has to use 
detectors provided with both polarizers and frequency 
filters. However, sometimes it is more convenient and, 
maybe, even more interesting to use only polarizers and 
wide-band detectors, non-selective in frequencies. In the- 
oretical description, such situation corresponds to aver- 
aging of the ququart's states over photon frequencies, or 
taking traces of the biphoton density matrix with respect 
to frequency variables of photons. In a general case, this 
gives rise to two-qubit biphoton mixed polarization states 
(MPS) 3]. In this paper features of such mixed states 
are briefly summarized and schemes for measuring their 
parameters are described. The method to be described is 
based on a series of coincidence measurements. In com- 
parison with the earlier suggested general scheme of mea- 
suring the ququarts's parameter (Appendix B of Ref. 0) 
the case of MPS has its own rather interesting peculiari- 
ties. 
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obey- 
ing the normalization condition = 1; 

a Ir h > a If i ' a v h i ana - a v i are one-photon creation oper- 
ators for four one-photon modes {H, u>h}, {H,u>i}, 
{V,0Jh}i an d {V,0Ji}. Superpositions of one-photon states 



where C< are arbitrary complex constants 
the normalization condition ^ i \Ci\ 2 



a' H h \0), a' H JO), Oy h |0), anddyJO) form one-photon qu- 
dits with the dimensionality of the one-photon Hilbert 
space d = 4. Double population of different one-photon 
modes corresponds to the basis state vectors in Eq. (QJ, 
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JO), etc. These basis state vectors, as well as 



their superpositions, can be considered as describing two- 
qudit states belonging to the two-photon Hilbert space 
of the dimensionality D = d 2 = 16. 

It's very fruitful to use not only state vectors of bipho- 
ton ququarts but also their wave functions. In a general 
case of arbitrary bipartite states with arbitrary variables 
of two particles x\ and xi the bipartite wave function 
^(x\,X2) is defined via the bipartite state vector |W) as 
^(x\,x-2) — (x±, X2\^}- One reason why it's important 
to use wave functions is that in terms of wave functions 
one can use the simplest definition of entanglement . Ac- 
cording to this definition a bipartite state is entangled if 
its wave function cannot be presented in the form of a 
product of two single-particle wave functions 



*(xi,x 2 ) ^ ip{x 1 )x(x 2 ). 



(2) 



II. BIPHOTON QUQUARTS AND MIXED 
POLARIZATION STATES 



In a general form, the state vector of an arbitrary 
polarization-frequency biphoton ququart is given by a su- 



Otherwise, if one can find such function f(xi) and x( x 2) 
that ^(xi,X2) — <p(xi) x( x 2), the state is disentangled. 

In the case of biphoton states with two degrees of 
freedom for each photon, for defining two-photon wave 
functions one has to introduce two pairs of variables, 
x\ = {crijOJx} and x 2 = {02,^2} • As photons are indis- 
tinguishable, we cannot attribute variable numbers 1 and 
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2 to any of two photons, though we know for sure that 
the amount of variables equals the amount of photons 
(2) times the amount of degrees of freedom (2), which 
gives 4 or two pairs. The biphoton wave function corre- 
sponding to the state vector of Eq. ([J} can be found with 
the help of the general rules of quantum electrodynam- 
ics (see, e.g., [H). The result can be written in different 
forms. The form most convenient for the further con- 
sideration is that related to the use of polarization and 
frequency Bell states 

* (4) (cti,wi; cr 2 , w 2 ) = * (3) (cri,o- 2 )*+(wi,a;2) 

+B^_(a 1 ,a 2 )^-(u} 1 ,ui 2 ), (3) 

where \&' 3 )(oi, a 2 ) is the wave function of a purely polar- 
ization qutrit 

f (3) (CTi,o- 2 ) = Ci* ffH (0-1,02) + B + ^ + (a 1 ,a 2 ) 

+C^ vv (ai,a 2 ) (4) 
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^ HH (o"l, 02) = 5 ai ^ H 5 a . 2 ^ H 
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and the Bell-state wave functions depending on polariza- 
tion or frequency variables are given by 

Sa 1 ,HSor 2 ,V ± S<7 1 ,V$a 2 ,H 
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superscripts "poP and "/r" in Eqs. ©-© indicate po- 
larization and frequency degrees of freedom. Besides, as 
it's clear from comparison of the functional and matrix 
forms of the biphoton wave functions in Eqs. (IS])-©, the 
upper lines in two-line columns correspond to the hori- 
zontal polarization and higher frequency uih and the lower 
lines - to the vertical polarization and lower frequency u>i . 
The constants B± in Eqs. (J3j) , (j4]) are expressed via C 2 ^ 
of Eq. <jTJ) as 



B. 



C 2 ± C3 
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(9) 



with the normalization condition taking the form \C\\ 



Note that the ququart's wave function ([3]) contains 
both symmetric and antisymmetric wave functions of the 
polarization and frequency Bell states. But the antisym- 
metric Bell-state wave functions ^ p ° l and \f r _ appear 
only in the form of their product, which makes the total 
wave function Vl^ 4 ) symmetric with respect to the trans- 
position of photon variables 1^2. This is the obligatory 
feature of two-boson pure states, which often is not taken 
seriously but which manifests itself, e.g., in the existence 
of MPS discussed below and in Ref. 3|. Besides, ow- 
ing to the symmetry, both terms in the ququart's wave 
function ^ obey the entanglement criterion @ and, 
hence, all biphoton ququarts are entangled. In some 
cases this is a purely frequency entanglement (e.g., if 
C4 = B + = B- = and C\ = 1), but in a general 
case entanglement of biphoton ququarts is an inseparable 
mixture of the polarization and frequency entanglement. 

The density matrix the state (TTJ), ^ is given by 



p(*) = *W®$Wt. 



(10) 



This density matrix characterizes pure states. But be- 
ing averaged over frequency variables, turns into the 
density matrix of a mixed two-qubit polarization state 
0. Written down in the basis {^ H h, \ 
and with dropped 12 zero lines and columns , the aver- 
aged density matrix takes a rather simple form 



Tr o (4) = P 



(3) 





(\Ci\ 2 




dB* 

\B+f 
CiB* 




CxCt 
B+Cl 
|C 4 | 2 






\B-\ 2 




\B-\V 



(11) 



(12) 



where 



p (3) = $(3) g)^(3)t 



(13) 



|B + | 2 + |£U 2 + |C 4 | 2 = 1. 



is the qutrit's coherence matrix [5|. In a general case 
the density matrix p characterizes a mixed polarization 
state. The only two exceptions occur in the cases -B_ = 
and |£?_| = 1. In the first case the ququart is reduced 
to qutrit, and in the second case the qutrit's contribu- 
tion to the quqaurt's wave function © equals zero (as 
Ci = C 4 = B+ = 0). In both cases S_ = and \B-\ = 1 
the ququart's wave function \E'( 4 ) factorizes for parts de- 
pending on the polarization and frequency variables sep- 
arately. This is the reason why in these cases averaging of 
a pure polarization-frequency state over frequency vari- 
ables leaves the remaining polarization state pure. In 
all other cases (\B-\ ^ 1, 0) there is no factorization for 
frequency and polarization parts in and, hence, the 
state, arising after averaging over frequency variables, is 
mixed. 

There are other forms of presenting the averaged polar- 
ization density matrix p alternative to that of Eq. (|12l) . 
One of them used below consists in the presentation of p 
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in the form of a sum of products of 2 x 2 single-photon 
polarization matrices: 



\Ci\ 



1 




2 




1 



1 




\B + \ 2 + \B-\ 2 
2 



1 




\Ca 




1 




1 



As it's clear from the definition of p, all matrices in this 
equation and further below refer to the polarization de- 
gree of freedom, with averaging over the frequency vari- 
ables already performed. For this reason, to shorten 
notations, here and below we drop the superscript pol 
common for all arising matrices. In Eq. (| 14[) only four 
products of 2 x 2 matrices are shown explicitly. In these 
four products all matrices are diagonal, whereas in all 
other 12 products, at least one of the matrices (* *). or 
( * * ) 2 is off-diagonal. Such terms do not contribute to 
conditional probabilities analyzed below in section V. 

The density matrix of MPS can be further reduced over 
polarization variables of one of two photons to give rise 
to the mixed-state reduced density matrix of the form Q 



ClB + +B* + C'^ 
71 



CiB* + +B+Cl 

|C 4 | 2 



%/2 

|B+| 2 +|B-| 2 
2 



(15) 



III. 



CORRELATIONS IN MIXED BIPHOTON 
POLARIZATION STATES 



Two correlation parameters found in the general form 
from the density matrices (|12l) , (|14l) , (|15l) are the Schmidt 
parameter K and concurrence C [3j: 



K 



l + (l-|B-| 2 ) 2 -|2CiC 4 -B| 



2 12 



and 



C 



\2dC 4 - B , 



IB- 1 



(16) 



(17) 



The concurrence C (jTTJ) characterizes the degree of en- 
tanglement or the amount of quantum correlations in 
MPS. Another quantifier of quantum correlations in such 
states is the so called relative entropy @ defined as the 
"distance" between the density matrix p and the density 
matrix a of the closest disentangled state 



S re i = Tr[p(\og 2 p - log 2 a)] 



(18) 



For MPS with C\ = C4 = the relative entropy was 
found in Q and shown to be less than concurrence at 
any values of the remaining nonzero parameters |~B-| 
and |£> + | = 1 — \B-\. The only exceptions occur at 
\B-\ = 0, landl/v^, where the concurrence and rela- 
tive entropy are equal. Thus, it was found that S re i < C, 



which can be interpreted as indication that the relative 
entropy is a better entanglement quantifier than concur- 
rence and that the latter can exaggerate slightly the de- 
gree of entanglement in the case of mixed polarization 
states. For such states, in accordance with the ideas of 
Refs. 6-8] one can define the quantifier of classical cor- 
relations as the difference between the von Neumann mu- 
tual information / = 2S(j) r ) — S(p) and relative entropy 



(14) S re l 



Ccl — I — S; 



rel • 



(19) 



As for the Schmidt parameter of mixed states K 
(|16p . in contrast with pure biphoton polarization states 
(qutrits) , K is not related anymore to the concurrence C 



CEZ): C^J2ll-K 1 ) and K ^ 1/(1 - C 2 /2) as in 



the case of pure states of biphoton qutrits. On the other 
hand, the Schmidt parameter of MPS remains related to 
their degree of polarization P by the same relation as in 
the case of pure states of biphoton qutrits 



P 



2 (l-K 



1. 



(20) 



where P = \S\, S = Tr(ap r ) is the vector of Stokes 
parameters, and a is the vector of Pauli matrices. Evi- 
dently, Tr {ap r ) ee Tr (a p^) and, hence, P = P^\ i.e., 
the degree of polarization of the mixed state coincides 
with the degree of polarization of the original two-qudit 
ququart, and they both are determined by the Schmidt 
parameter of the mixed state K via Eq. ([2U|) . 

As the degree of polarization is a classical concept, we 
can deduce from Eq. (|2"0)) that in the case of mixed states 
the Schmidt parameter K is related to the amount of clas- 
sical rather than quantum correlations. In terms of K, 
a new parameter characterizing the amount of classical 
correlations can be defined as 



Cd = \ 2 l-K 



(21) 



It may be interesting to notice that for the state with 
C\ = Ca = this parameter coincides with that of (fT9]) 

C c i = Cd. (22) 

Note finally that in other special cases, £?_ = or 
\B- \ = 1, when states averaged over frequencies remain 
pure, all discussed parameters of quantum and classical 
correlations coincide with each other and are equal to a 
half of the von Neumann mutual information 



C — Sr 



C c l — C c 



1/2. 



In these cases the relation C = 



1 -K 



(23) 



becomes 



valid again, and this is the reason why in pure bipartite 
states the Schmidt parameter K can be used for charac- 
terization of the amounts of both quantum and classical 
correlations. 
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IV. COMPARISON WITH A TWO-QUBIT 
PURE-STATE MODEL OF BIPHOTON 
QUQUARTS 

A picture of mixed two-qubit polarization states de- 
scribed above differs significantly from traditionally 
widely used model of two-qubit pure-state ququarts. 
This model starts from the same state vector as given 
by Eq. ([J). But then frequencies of SPDC photons u>i 
and u>2 are considered as given numbers rather than vari- 
ables, e.g., as cji = ujh and cu2 = uji- Owing to this, 
two photons of SPDC pairs are considered as "partially 
distinguishable" , owing to which the polarization bipho- 
ton wave function appears to be not necessarily symmet- 
ric with respect to the transposition of particle's variables 
1^2, and can be written in the form 

^feC "!.^) = *cn? (0-1,02) +B + y + (a 1 ,a 2 ) 

+C i ^ vv (a 1 ,a 2 ) + B-^-(ai,a 2 ). (24) 

This is a wave function of a pure two-qubit state, and it 
yields the well known results for the Schmidt parameter 
and concurrence: 

Kf qh = (25) 

G% = ^2 (l-K~} b ) =\2CiC A -B 2 + +B 2 _\. (26) 

In a general case, these expressions differ from K and C 
of Eqs. (fl"6)l and flTTl) . We believe that the correct results 
are those based on the picture of two-qudit polarization- 
frequency bipohoton ququarts and of MPS arising after 
averaging over frequencies, i.e., the results determined by 
Eqs. (fT6|) -pO ]) . Weak points of the two-qubit theory of 
biphoton ququarts are evident. Photons of SPDC pairs 
are always indistinguishable. If there is something that 
looks like "partial distinguishability" , this is an indica- 
tion that there is, in fact, an additional degree of freedom, 
and with this degree of freedom taken into account ac- 
curately, photons are evidently indistinguishable. Wave 
functions of two photons in a pure state cannot be asym- 
metric with respect to the transposition of their variables. 
Its symmetry is dictated by the Bose-Einstein statistics 
of photons. This feature is clearly violated in the two- 
qubit wave function of Eq. (|24[) where the symmetric and 
antisymmetric Bell-state wave functions are summed on 
equal terms. Note however, that a simple symmetriza- 
tion of the expression in Eq. (f2~4"| vanishes the term, pro- 
portional to VP-, and reduces the ququart's wave func- 
tion *2gi t0 tliat of a qutrit * (3) g]). To get a correct 
quqaurt's wave function ([3]), in addition to symmetriza- 
tion, one has to give freedom to photon frequencies wig 
by considering them as variables which can take one of 
two values each: either cui = lo^ and u-i = U)i or cjj = u>i 
and 102 = ^>h- Actually, this means that we never know 
which photon has which frequency. Averaging of states 
of bophoton ququarts over frequencies gives rise to 



MPS considered here and in Ref. Q. There is no way to 
get such states in a two-qubit model. In principle, differ- 
ences between predictions of the theory of mixed states 
and of the two-qubit model of ququarts can be seen in 
experiments on measurement of the degree of polariza- 
tion of biphoton polarization-frequency ququarts. Some 
simple examples of experimental schemes where these dif- 
ferences are well pronounced are described in Ref. [3|. 



V. RECONSTRUCTION OF QUQUART'S 
PARAMETERS IN EXPERIMENTS 

The next questions are how to measure in experiments 
parameters of MPS and of pure states of polarization- 
frequency ququarts. It was shown earlier Q that, in 
principle, a series of coincidence polarization-frequency 
measurements in three different bases provides sufficient 
amount of data to get a complete set of equations for find- 
ing all ququart's parameters. But equations obtained in 
such a way were rather complicated and not convenient 
for practical purposes. Here we will consider first a sim- 
pler problem of finding parameters of the above discussed 
MPS related to biphoton ququarts. And then, at the 
last stage, we will show how this procedure can be pro- 
longed in a very simple way to reconstruct explicitly all 
ququart's parameters. Note also that the methods of us- 
ing series of coincidence measurements for reconstructing 
parameters of quantum states are alternative to standard 
and rather widely used methods of quantum tomography 
for biphoton ququarts [2, • 



A. Independent constants characterizing biphoton 
qutrits, ququarts and mixed polarization states 

Pure states of qutrits (Q| and ququarts © are char- 
acterized, correspondingly, by three and four complex 
parameters, {C\, B + , C4} and {C\, B + , C4, -B-}, which 
corresponds to 6 and 8 real constants. But these param- 
eters are not completely independent: there are normal- 
ization conditions and, besides, in both cases the common 
phases of wave functions do not affect measurable quan- 
tities and, hence, can be taken having arbitrary most 
conveniently chosen given values. These conditions re- 
duce the amount independent real constant parameters 
characterizing qutrits and ququarts, correspondingly, to 
4 and 6. MPS considered above occupy an intermediate 
position between qutrits and ququarts. Parameters char- 
acterizing these states are the same as in the case of pure- 
state ququarts, {C\, B + , C4, B-}. But, in addition to 
the normalization and common phase conditions we have 
now one condition more: as seen well from the structure 
of the density matrix p written in the form (1121) , features 
of mixed states do not depend of the phase of the parame- 
ter B_ . This is seen well from the derived expressions for 
the Schmidt parameter and concurrence (1161) and (|17l) . 
which depend on _B_ only as on Also, as seen well 
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from the definition of p in the form of Eqs. (ITT1) . (|T3)) . the 
density matrix of MPS does not depend on the phase of 
the qutrit's wave function which enters into the def- 
inition of p as determining one of its components. Owing 
to this, the phase of can be chosen, e.g., in a way, 
making the parameter B + real and positive. Thus, we 
find that in this case MPS are characterized completely 
by 5 independent real constants: \C\\, tp%, B + , | C4 1 , 994 , 
where ^1,4 are phases of the parameters 6*1,4 with the 
constant \B—\ to be found from the normalization condi- 
tion. 



B. Conditional probabilities and coincidence 
measurements 

By definition, the conditional probability Wo-I^z is the 
probability for a photon 1 to have polarization a under 
the condition that the second photon (2) of the same pair 
has polarization a' . Relations between the conditional 
probability and parameters of MPS follow directly from 
the presentation of the density matrix p in the form of a 
sum of products of 2 x 2 single photon matrices (fT4|) : 



w H \ v = w v \ H 



W V \y = 



|C 4 | 2 , 



(27) 
(28) 



Owing to normalization, ^2 a cr , w a \ a , = 1. 

In experiment, conditional probabilities can be found 
from coincidence measurements. For this goal a bipho- 
ton beam has to be divided for two channels 1 and 2 by 
a non-selective Beam Splitter (BS). Half of SPDC pairs 
will be divided between channels, whereas another half 
of undivided pairs will appear either in the channel 1 or 
2 (see a scheme in Fig. [T]). In coincidence measurements 



V« A. y(V) 





x(H) 



Figure 1: (a): a scheme of experiment for coincidence measure- 
ments, BS - beam splitter, P - polarizers, D - detectors, M - mirror; 
(b): the horizontal- vertical coordinate frame and the frame turned 
of an angle a. 

only divided pairs are registered, and for such pairs num- 
bers of channels 1 and 2 can be associated with photon 
or variable numbers 1 and 2 in Eqs. ©-© an( i (O- 
Measurements consist in counting photons with a given 
polarization. Selection of polarization is provided by po- 
larizers installed in each channel in front of detectors. 
Orientation of each polarizer can be changed indepen- 
dently from horizontal to vertical and vice versa. Each 
series of measurements with given polarizer orientations 



has to be performed under identical conditions and to 
take the same time. As the undivided pairs do not partic- 
ipate in such measurements, the amount of photons with 
any given polarization accessible for registration is twice 
less than in the original beam. Besides detectors have 
some efficiency less than 100%, which further diminishes 
amounts of registered photons. But for relative amounts 
of registered photons all these losses do not matter. A 
computer obtaining signals from both detectors registers 
only coinciding incoming signals and provides measure- 
ments of the relative amounts of counts coinciding with 
conditional probabilities ([2"7). (j2"gj) . Let N a \ a , be the 
counted number of pairs with the photon polarization a 
in the channel 1 and a' in the channel 2. Then we have 



Z-/(T,< 



N„ 



(29) 



Now Eqs. ([27]) and ([29]) can be used for finding two real 
constants characterizing MPS, |Ci| 2 and | C 4 1 2 , directly 
from experimental measurements with polarizer in both 
channels 1 and 2 oriented either horizontally or vertically. 
Measurements with other orientations (HV) and (VH) 
are necessary too but only for determination of the nor- 
malizing factor in the denominator of Eq. (|29p . The sum 
\B+\ 2 + |B_| 2 = 1 - |d | 2 - |C 4 | 2 is determined by nor- 
malization but the constants B+ and |£?-| 2 themselves 
remain undefined, as well as the phases tpi and tp^. For 
finding them other measurements are needed. 



C. Conditional probabilities in rotated frames 

Additional information about parameters of MPS can 
be obtained from coincidence measurements with polar- 
izers in both channels 1 and 2 turned for the same angle a 
with respect the horizontal or vertical axes. The results 
of such measurements (photon counting) are related to 
the corresponding conditional probabilities by the same 
relation as in the case a — (|29p : 



N a \ a , _ , N 90 . +a \ a , _ 
- w a \ a , , , — w 90 ° +a \ a , 

90° and N tot 



N, 



(30) 



tot 

where a' = a or a + 90° and N to t = J2 a '{ N <x\a' + 

N9Q°+a\ a i] 

For finding the conditional probabilities w a \ a , and 
wgo°+ Q | Q /, we have to rewrite the same ququart's wave 
function as given by Eq. ^ in the frame turned for an 
angle a around the photon propagation axis Oz (in the 
(x a , ?/ Q )-plane in Fig. 16). Transformation to this frame 
is provided by the basic transformation formulas for one- 
photon states 



cos a 







^ — sin a 












) + cos a 









(31) 
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where (J)" and (°) correspond to polarizations along 
directions a and 90° + a. Evidently, the transformations 
(|3Tj) do not affect the frequency part of the ququart's 
wave function ([3]). Besides, as can be easily checked, 
the polarization antisymmetric Bell-state wave function 
is invariant with respect to the transformations pip , i.e., 
'J" expressed via (p) and (°) has the same form as 

ty- expressed via (q) and (°) [Eq. This means that 
after the transformation (|3ip the part of the ququart's 
wave function with the product of antisymmetric Bell 
states does not mix up with the symmetric part, and the 
form of ququart's wave function is invariant with respect 
to transformation: 



$( 4 ) = $( 3 ) Q ^{ r + B^ p _° la ^ f _ r 



(32) 



where \I/( 3 ) a has the same form in the transformed basis 
(l) a and (°) a as * (3) in the original HV basis (gj) 



^(3) a = Qa 9 +B «^* + d * 90 o 



+a,90°+a- 



(33) 



The coefficients Cf , B°l_ and C" are easily found to be 
given by 2] [Eq. (A. 13)] 

C? = cos 2 a cos a sin a -B_|_ 4- sin 2 (1C4, (34) 

B% = -V2 cos a sin a (C\ — C 4 ) + cos2a-B + , (35) 
C 4 = sin 2 a C\ — V% cos a sin a £> + + cos 2 a C 4 , (36) 

The conditional probabilities in the turned frame are de- 
fined similarly to their definition in the HV- frame (|27[) : 

v> a \ a = \C?\ 2 

= |cos 2 aCi + v / 2cosasinai? + +sin 2 aC 4 | 2 , (37) 

W90°+a| 90 o +Q = IC" | 2 

= I sin 2 aCi - v2 cos a sin a B + + cos 2 a C4 1 2 , (38) 

W <*\90°+ a = w 90°+a| Q = Z • (39) 



If a is small, in the linear approximation in a, Eqs. ([371 
and (1381) take the form 



w a \ a = w H \ H + 2V2 a\Ci\B+ cos ip u (40) 

ti>90°+a| 90 o +Q = w v \ v - 2v / 2a|C , 4|S + cos^ 4 , (41) 

where, as assumed, B + is taken real and positive. By 
denoting 



lim^o w ^~ a WH ^ =tang 1 , 
lim^o + ' 90 ° + ° — = tan 04 



we can rewrite Eqs. (|40l) and (14"TT) as 



cos </?i 



tan $i 



2 3 / 2 |C 1 |B, 



COS <^4 



tan 9a 



2 3 / 2 |C 4 |S+' 



(42) 



(43) 



Here expressions for the parameters \C\\ and | C 4 1 in 
terms of conditional probabilities are known (|27[) . and 



the derived Eqs. (|40l) and (|4T|) determine directly phases 
ipi and if4 as functions of For finding B + [and then 
from Eq. ([28]) ] one has to make one measurement 
and one derivation more, for example, in the frame with 
a = 45°, i.e., in the frame turned for 45° with respect to 
the original HV- (xy-) frame. In this case we find from 
Eq. |J35J) that Bf° = \C X - C 4 | 2 /2 and, hence, Eqs. 
and (|39|) yield the following equation: 



w 45 o | 1350 



\C X - C 4 



B\ 



wh\ v --±, (44) 



where 



\d-dl 2 = |C 1 | 2 + |C 4 | 2 -2|C 1 ||C 4 |cos(^ 1 -( /34 ) (45) 
and 

tan 6\ tan 64 



cos(y>i - 934) 



2 3 / 2 |Ci|i? + 2 3 / 2 |C 4 |R, 



tan 81 \ 2 ( tan 64 



2 3 / 2 \d\B + J V2 3/2 |C 4 |fi 



1/2 



(46) 



With known Ci j4 and tan 951,4, the only unknown pa- 
rameter in Eq. (|44| [combined with Eqs. (|45|) and 
([4*61 ] is B + , and this equation has to be solved numer- 
ically. When B + is found, Eq. $ZE§ yields \B-\ = 
1/2 

[2 wh \y — l-B+l 2 ] ■ This concludes determination of 
all five parameters of MPS in terms of conditional prob- 
abilities. 



D. Scenarios for experimental measurement of the 
parameters of mixed polarization states 

In all cases, the first step consists in performing at 
least three series of coincidence measurements of photon 
numbers with polarizers in channels 1 and 2 installed 
along either horizontal or vertical directions, Nh\ h , 
Nh\ v , and Ny\ v (because of photon indistinguishabil- 
ity N V \ H = N H \ V ). Then, with the help of Eqs. (JSTJ) 
we find two constants, |Ci| and |C 4 | plus the relation be- 
tween l-B+l 2 and |-B-| 2 (l2"5j) . The next steps are different 
for the situations of zero or non-zero obtained values of 
the parameters \C±\ and |C 4 |. 



1. c 1 = c i = o 

If the above-described measurements with horizontal- 
vertical orientations of polarizers give |Ci| = |C 4 | = 0, 
the remaining two constants to be found are \B + \ and 
\B—\, and their measurement is very simple. E.g., from 
of Eq. (J211) we find Cf = sm2aB + /y/2. Then Eq. ([37]) 
yields 



IB+I 2 = 



2 w n 



2a 



2N a \ a /N tot 
sin 2 2a 



(47) 
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and |-B_| = y 1 — l-B+l 2 - Thus, if a complete set of 
horizontal- vertical coincidence measurements gives Ci = 
C4 = 0, for finding \B+\ and |£?_| one has to make only 
one coincidence measurement more, with identically ori- 
ented polarizers P\ and Pi and with arbitrary chosen 
angle of their orientation a 0, n. 



2. G x 0, C 4 = 

If the horizontal-vertical coincidence measurements 
give Ci 7^ 0, C4 = 0, at a = 45° the formula in Eq. 
(1351) is reduced to Bjf 5 = Ci/y/2, owing to which Eq. 
(l39ll takes the form 



iVl 35 o| 45 o _ ICil 2 , |B_| 2 



JVt, 



(48) 



This is the equation for finding \B_\, after which 
\B + \ is also easily found from normalization |B_| = 
a/1 — I Ci 1 2 — l-B+l 2 . Thus, for measuring absolute val- 
ues of all three constants, |d|, \B+\, and \B-\ ( with 
C4 = 0), it is sufficient to complete horizontal- vertical 
coincidence measurements by measurement of the coin- 
cidence number of photons with polarizers Pi and P2 
turned, correspondingly, for 135° and 45° with respect 
to the horizontal direction. In principle, in addition to 
these three constants there is one constant more charac- 
terizing MPS with C4 = 0, the phase ipi of the parameter 
C\ ( with real B + ). The way of its measurement is de- 
scribed in the following subsubsection. But it should be 
noted that in the case C4 = the MPS correlation pa- 
rameters K (fl~6|) and C (jTTJ) do not depend on tp\. Note 
also that the described scheme of measurements is valid 
also in the case C\ = 0, C4 7^ 0. 



3. Nonzero Ci,4 

Let now both constants \C\\ and IC4 1 found from the 
horizontal-vertical coincidence measurements be differ- 
ent from zero, C1.4 7^ 0. Then, the procedure of mea- 
suring other MPS parameters is more complicated. In 
particular, for measuring phases ipi^ of Ci^, we suggest 
to use their relations (|4"3"|) with the parameters tan 6*1,4 
(|4"2")l characterizing the behavior of the function w a \ a (a) 
in small vicinities of the points a = and a = 90°. 
Specifically, for finding B\ and <pi we suggest to measure 
coincidence numbers of photons N± ao \ ±ag with both po- 
larizers Pi and Pi turned for some small angles olq and 
— olq with respect to the horizontal direction (e. g., with 
Q'o = 5° = 0.087 rad -C 1). In accordance with Eqs. 
pop and together with the earlier made measurement of 
wh\h tms gives three values of the function w a \ (a) at 
three values of a, a = —ao, 0, andao- In Fig. [5] these 
values are indicated by letters A, O, and B. With this 
three values we can reconstruct the function w a \ a {a) 
in a small vicinity of the point a = 0. The pictures 




tOfl 1 u O 


A \ 












1 a 



(c) A w a \ 



w H \ H 


V 


A 













1 a 



-ao a a 

Figure 2: Conditional probability w a \ a (a) at \a\ <C 1 (solid lines) 
and results of its measurements A, O, and B. 



(a) , (6), and (c) correspond to three different possible lo- 
cations of the points A and B. The picture (a) corre- 
sponds to the case when the points A and B are sym- 
metric with respect to O, and all three points A, O, and 
B can be connected by a single straight line. In this case 
the angle ao is small enough for validity of the approxi- 
mation linear in a in Eq. f|40[) in all range [— cto,ao\. In 
the case (b) positions of the points A and B are asymmet- 
ric, and the points A, O, and B can be connected only by 
a parabola. In this case, the straight line corresponding 
to the linear approximation of Eq. (|4T)|) is tangent to the 
parabola in the point O. The angle 9\ is determined in 
these two cases as the angle between the horizontal line 
and either the line w a \ (a) in the case (a) or the line 
tangent to the curve w a \ a (a) in the point O in the case 

(b) . In both cases, with known \C\\ and 81, we can use 
Eq. (|43l) for finding the phase (pi of the parameter C\ as 
a function of B + . Similar measurements and calculations 
can be done with polarizers deviating for a small angle a 
from the vertical orientation, to determine #4 and, then, 
ip4 as a function of B + . The last step is a single coin- 
cidence measurement with one of two polarizers turned 
for 45° and the other one for 135°. This measurement 
gives wi35o| 450 , and then Eqs. P3 )) -([lB ]l can be used for 
finding B + and, then, tpn- With \B—\ found from Eq. 
(|28[) . this finalizes the procedure of finding all parameters 
of MPS of a general form. 



4. B+=0 

In the picture (c) of Fig. [5] the points A and B are 
located symmetrically with respect to the vertical axis 
crossing the point O. Again, three points A, O, and B 
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can be connected by a parabola, 

w a \ a w w H \ H + ka 2 . (49) 

The line tangent to this parabola in the point O is hori- 
zontal. This means that in this case there is no validity 
region for the linear approximation of Eqs. (j40|) and fill) , 
which is possible (at 7^ 0) only if B+ = 0. In this 
case a direct measurement of the parameter k in Eq. (|49[) 
appears to be sufficient for completing the reconstruction 
of the parameters characterizing MPS. Indeed, in the case 
B + = and a < 1 Eq. (gD]) takes the form 

C? « Ci + (C 4 - Ci) a 2 . (50) 

Moreover, as B + = 0, wc can choose now an arbitrary 
phase of in a way providing <pi = 0, which leaves 
only one phase parameter </? 4 to be determined form ex- 
periments in addition to C\ and C4. Under this assump- 
tion Eq. (f50|) gives 

|C?| 2 « Idl 2 -!- 2|C 1 |(|C 4 |cos ¥ > 4 -\C X \) a 2 (51) 

and, being compared with Eq. (|49|) . 

fc = 2|C*i|(|C* 4 |cos^4-|Ci|). (52) 

Thus, by measuring experimentally the parabola param- 
eter k of Fig. 2(c) and Eq. (JSDJ we find from Eq. ([52]) 
the phase ^4, and this concludes the reconstruction of all 
parameters of MPS in the case £?+ = 0. 

VI. RECONSTRUCTION OF THE QUQUART'S 
PARAMETERS 



By assuming again that tp + — 0, we find from this equa- 
tion cos <p_ expressed in terms of the experimentally mea- 
surable relative amounts of photon counts. 

The case B + = (but Ci,4 7^ 0) is not much more 
complicated or difficult. With polarizers in both chan- 
nels 1 and 2 turned for an arbitrary but identical an- 
gle a, one can use, in fact, the same scheme of mea- 
surements as described above for the case B + = 0. In- 
deed, as mentioned above, in the case B + = Eq. (|35|) 
takes the form B" = v2 cos a sin a (C 4 — C\). As the 
parameters C\ (with pi = 0), and C4 = \C^e lLpi are 
supposed to be known already from purely polarization 
measurements, we can write the difference C4 — C\ as 
C 4 - Ci = |C 4 - Cile*'^*- 1 , where |C 4 - Ci| and 
are easily calculable. Now, with the turned polarizers, 
one can measure the conditional probability of register- 
ing a high-frequency photon polarized in the direction a 
in the channel 1 and a low-frequency photon polarized 
in the direction a + 90° in the channel 2. This condi- 
tional probability is related to the unknown phase (/?_ by 
a formula very similar to that of Eq. ((55|) 

w *,h\ a+9 o°,i - Ntot ~ 2 

\B«\ 2 + |5_| 2 + 2\B°\\B_\ cos(^_ - ^4-1) 
= — ^ , (54) 

where \B%\ = -^| sin2a||C 4 - C x \. Eq. fjlj) can be 
used for finding the phase ip- from the data to be ob- 
tained from the coincidence polarization-frequency mea- 
surements in the case of ququarts with B+ = 0. 



Let us return now to pure states of biphoton ququarts. 
The analysis of the previous Section shows that by means 
of purely polarization measurements one can determine 
all parameters of the "qutrit's part" of the ququart 
(except its phase) and the parameter \B_\. The only re- 
maining unknown parameter of the ququart 's states is 
the phase p- of the parameter _B_ . Though characteris- 
tics of MPS do not depend of y>_, features of pure states 
of ququarts can be phase-sensitive. The phase ip_ cannot 
be found from any purely polarization measurements and 
requires combined polarization-frequency measurements. 
This means that in the experimental scheme of Fig. [2] 
one has to install in front of detectors both polarizers 
and frequency filters. Such coincidence measurements 
sufficient for determining the phase <p- are most simple 
in the case of ququarts with B + 7^ 0. Then, one of the 
measurable conditional probabilities is the probability of 
registering high-frequency horizontally polarized photons 
in the channel 1 under the condition that simultaneously 
one registers low-frequency vertically polarized photons 
in the channel 2 

N H ,h\ V ,l \B++B-\ 2 

WH ' hlv ' i= ^r = — 2 — 

_ \B+\ 2 + |B_| 2 + 2\B+\\B_\ cos(p- - p+) 



VII. CONCLUSION 

Thus, biphoton ququarts are more complicated and 
their physics is more rich and interesting than usually 
assumed. The key elements of this newer understand- 
ing are (i) the obligatory symmetry of biphoton wave 
functions (in pure states) as wave functions of two indis- 
tinguishable bosons, and (ii) consideration of frequencies 
of photons in biphoton polarization-frequency ququarts 
as variables independent of polarizations rather than as 
given numbers. In this approach biphoton polarization- 
frequency ququarts are states having two degrees of free- 
dom for each photon, polarization and frequency Owing 
to this, all biphoton polarization-frequency ququarts are 
entangled and, in a general case, their entanglement is 
an inseparable mixture of the polarization and frequency 
entanglement. Another interesting consequence of the 
formulated features of biphoton ququarts concerns their 
images to be seen in experiments. If in fully polarization- 
frequency coincidence measurements ququarts are seen as 
pure states, in simpler purely polarization (non-selective 
in frequencies) measurements the same states are seen as 
two-qubit mixed polarization states, MPS. MPS are char- 
acterized by the ququart 's density matrix reduced with 
respect to the frequency variables. Parameters of MPS 
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are found. They appear are to be rather peculiar, dif- 
fering significantly from those of the full-dimensionality 
ququarts, and rather useful. In particular, the Schmidt 
parameter of MPS, K is found to be related directly to 
the degree of polarization of ququarts. Features of MPS 
can be used straightforwardly for experimental measure- 
ment of ququart's parameters if they are not known in 
advance. A scheme of such measurements is suggested 
and described. The main idea is in separation of a bipho- 
ton beam for two channels by a simple non-selective beam 
splitter and in performing series of coincidence measure- 
ments with photon counters. The first stage consists in 
making purely polarization measurements with different 
orientations of polarizers in both channels and in find- 
ing in this way all parameters of MPS. Then for recon- 
structing completely the ququart's state, one has to find 
additionally only one of its phase parameters, which 



requires making one additional polarization-frequency 
coincidence measurement with both polarizers and fre- 
quency filters installed in both channels in front of de- 
tectors. So, the scheme suggested for reconstruction of 
the ququart's parameters separates purely polarization 
and polarization-frequency measurements and minimizes 
the amount of polarization-frequency measurements. We 
believe that such experiments are feasible and their re- 
sults may be sufficiently interesting. 
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